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Analysis  of  Air  Breathing  Hall-Effect  Thruster 


Leonid  Pekker1 

ERCJnc.,  Edwards  AFB,  CA  93524 
and 

Michael  Keidar2 

The  George  Washington  University,  Northwest  Washington,  DC  20052 


The  principle  idea  of  using  air  breathing  electrical  propulsion  for  a  vehicle  flying  at 
orbital  speed  on  the  edge  of  Earth’s  atmosphere  is  examined.  In  this  paper,  we  present  a 
simple  model  of  a  Hall  Effect  thruster  in  which  the  propellant  is  an  ambient  air.  The 
required  lengths  of  the  thruster  chamber,  the  magnetic  fields,  the  thrust,  and  other 
parameters  of  an  ideal  air  breathing  Hall  Effect  thruster  are  calculated  as  a  function  of  the 
flying  altitude  of  the  vehicle.  We  demonstrated  that  an  air  breathing  Hall  thruster  is  indeed 
capable  of  providing  the  needed  thrust  for  near  space  satellites  in  the  case  of  the  altitudes  of 
about  80  -  90  km. 


Nomenclature 

=  mass  of  a  heavy  particle,  in  au 
=  strength  of  the  magnetic  field 
Ddiff  =  ambient  gas  diffusion  coefficient 

E  =  strength  of  the  electric  field 

*  jet  =  energy  flow  leaving  the  Hall  thruster  with  the  plasma  jet 

r-  IOSS 

11 anode  =  heat  losses  at  the  anode 

jp  loss 

c-waii  =  energy  losses  at  the  wall 

Ejnnti  —  ionization  cost 

e  =  electron  charge 

Fjrag  =  drag  force 

fdras  =  drag  factor 

je  =  electron  current  density 

ji  =  ion  current  density 

H  =  altitude  of  satellite  orbit 

L  =  length  of  the  Hall  thruster  chamber 

In  A  =  Coulomb  logarithm 

M  =  mass  of  a  heavy  particle 

=  ion  mass 

me  =  mass  of  an  electron 

ne  =  electron  number  density 

nKas  =  number  density  of  ambient  gas 


1  Senior  Research  Staff  Engineer,  leonid.pekker.ctr@edwards.af.mil. 

2  Associate  Professor,  keidar@gwu.edu. 
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ion  number  density 
number  density  of  plasma 

electric  power  put  into  the  discharge 

inner  radius  of  the  opening  of  the  air  breathing  Hall  thruster 

outer  radius  of  the  opening  of  the  air  breathing  Hall  thruster 

Debye  radius 

electron  cyclotron  radius 

effective  cross-section  of  the  satellite  due  to  drag  force 

electron  temperature 

electron  temperature  in  the  plasma  jet 

ambient  gas  temperature 

ion  temperature  in  the  plasma  jet 

thrust 

Bohm  velocity 
electron  drift  velocity 
electron  thermal  velocity 
thermal  velocity  of  ambient  gas 

the  first  cosmic  velocity,  7.8  km/s 
coefficient  of  secondary  electron  emission 
delta-V  (change  in  the  velocity) 
ambient  gas  mean  free  path 
electron  mobility 

collision  cross-section  of  ambient  gas  molecules 
ionization  cross-section 
electron  lifetime  in  Hall  thruster 
ionization  time 

characteristic  time  for  gas  molecule  collisions  with  wall 

time  that  a  gas  molecule  is  flying  through  the  thruster  chamber 

electron-ion  collision  frequency 

the  Bohm  sheath  potential 

applied  voltage 

electron  cyclotron  frequency 


I.  Introduction 

This  work  is  motivated  by  increasing  interest  in  military  and  civil  spacecrafts  flying  at  the  altitudes  of  70  -  120 
km.  Since  the  drag  at  such  altitudes  is  still  significant,  the  thrusters  of  such  a  satellite  should  work  continuously  in 
order  to  maintain  the  orbit  altitude.  This  demands  that  a  significant  amount  of  propellant  be  stored  on-board 
satellites  that  use  ordinary  thrusters  for  propulsion.  However,  for  such  spacecrafts,  using  air  breathing  thrusters,  in 
which  the  propellant  is  ambient  air,  looks  very  attractive,  particularly  for  long-duration  missions;  air  breathing 
thrusters  allow  significantly  increased  payload-to-weight  ratio  for  such  satellites.  Using  air  breathing  thrusters  in 
Mars  and  other  planets  atmosphere  becomes  even  more  attractive  because  of  multi-years  satellite  missions  and  high 
cost  of  payloads. 

In  a  recent  paper  [  1],  Diamant  has  proposed  a  two-stage  cylindrical  Hall-Effect  thruster  for  air  breathing  electric 
propulsion,  in  which  the  first  stage  is  an  electron  cyclotron  resonance  ionization  stage  and  the  second  stage  is  a 
cylindrical  Hall  thruster.  The  author  built  such  a  thruster  and  demonstrated  its  operation  in  xenon  gas.  As  follows 
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from  his  assumptions,  his  two-stage  cylindrical  Hall  thruster  is  able  to  work  at  a  220  km  orbit  with  ambient  air 
passively  compressed  by  a  factor  of  500  [1].  However,  achieving  such  a  large  passive  gas  compression  seems  to  be 
difficult  and  will  be  associated  with  substantial  drag.  In  the  BUSEC  conception  of  an  air  breathing  Hall  thruster  [2], 
the  compression  of  air  is  achieved  by  a  diffuser.  It  is  not  clear  what  compression  ratio  of  the  ambient  air  can  be 
achieved  by  using  such  a  diffuser  at  orbits  of  about  100  km,  where  the  gas  mean  free  paths  are  tens  of  centimeters 
and  comparable  with  the  dimensions  of  the  thruster.  Furthermore,  the  use  of  a  diffuser  can  significantly  increase  the 
drag  because  gas  is  reflected  diffusely  in  space  [3]. 
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Figure  1.  Schematic  of  air  breathing  Hall  thruster  (not  to  scale). 


In  the  present  work  we  propose  a  simple  configuration  of  an  air  breathing  Hall-Effect  thruster,  (shown  in  Figure 
1),  in  which  the  incoming  air  flow  is  ionized  and  accelerated  directly  in  the  Hall  thruster  chamber  without 
preliminary  compression  as  in  [1,2].  Such  a  design  leads  to  a  decrease  of  the  possible  drag,  however  its  feasibility 
warrants  careful  consideration,  which  is  the  subject  of  this  paper.  The  description  of  the  model  and  numerical  results 
are  provided  in  Sections  II  and  111.  respectively,  and  the  conclusions  are  given  in  Section  IV. 

IT.  Description  of  the  Model 

The  following  assumptions  are  made  in  the  model:  (1)  no  drag  due  to  ambient  gas  passing  through  the  Hall 
thruster  chamber,  illustrated  in  Figure  1,  corresponding  to  a  condition  in  which  the  ambient  gas  freely  flows  through 
the  thruster  chamber  without  a  significant  interactions  with  the  wall;  (2)  the  gas  jet  leaving  the  thruster  chamber  is 
fully  ionized  providing  the  maximum  thrust;  (3)  the  plasma  in  the  Hall  thruster  is  quasineutral,  «,  =  n(.  =  np, , 

corresponding  to  a  condition  in  which  rD  is  smaller  than  the  characteristic  dimensions  of  the  thruster;  in  our  case. 
Figure  1,  rD  must  be  smaller  than  rh~ra ;  (4)  the  electron  cyclotron  radius  is  much  smaller  than  the  Hall  thruster 
gap,  rL  <  rb  -  ra  ,  as  is  required  for  confinement  of  electrons  in  the  Hall  thruster. 

The  first  assumption  is  well  satisfied  when  the  characteristic  time  required  for  a  gas  molecule  to  reach  the  wall 
and  transfer  its  momentum.  TwaU  ,  is  larger  than  the  time  required  for  the  molecule  to  freely  fly  through  the  thruster 
chamber,  rCl  =  L! V0 ,  i.e. 


Twull  •  ^1) 


L 


>1  . 


(1) 


In  the  case  of  free  molecular  flow,  when  the  gas  mean  free  path  is  larger  than  rb-ru,  Figure  1 ,  Twall  can  be 
estimated  as 


r  rb~ra 

Lwall  . . 

VTgas 


(2) 
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and  in  the  case  of  coliisional  molecular  flow,  rb-ra<  Agas ,  as 


Tvt -all  - 


Here,  we  have  substituted  ( ngas 


(rb~ra)~  (rb~ra)~  ._ngas'<7ijas'(rb  ra ) 

^diff  ^gm  ^Tgas  ^Tgas 

■  a  for  Agas  .  Combining  Eqs.  (2)  and  (3)  we  obtain 


(3) 


*wu  =  Max 
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The  second  assumption  is  fulfilled  when  the  ionization  time  of  a  gas  molecule, 
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is  much  smaller  than  r0 .  Let  us  estimate  Tioniz .  The  equation  describing  the  mass  conservation  law  for  the  heavy 
particles  can  be  written  as 


■  ngas  ripi 


Vi  + 


2  -ftp 

M, 


(6) 


where  the  left-hand  side  in  Eq.  (6)  is  the  flux  of  incoming  gas  and  the  right-hand  side  is  the  ion  flux  leaving  the 
chamber.  In  Eq.  (6)  we  have  assumed  that  the  downstream  plasma  in  the  thruster  chamber  is  quasineutral  and  fully 
ionized. 

It  should  be  stressed  that  in  the  model  we  assume  the  ambient  gas  is  pure  molecular  nitrogen.  This  is  a 
reasonable  assumption  since  the  air  contains  78  percent  of  nitrogen  and  21  percent  of  oxygen  and  their  electron 
impact  ionization  cross  sections  are  similar  [4].  In  the  model,  we  neglect  the  dissociation  of  molecular  nitrogen 
assuming  that  M  =  M , .  Since  for  electron  energies  between  20  -  80  eV  the  ionization  and  dissociation  cross 
sections  of  molecular  nitrogen  are  close  and  the  ionization  cross  section  of  molecular  nitrogen  is  much  larger  than 
the  ionization  cross  section  of  atomic  nitrogen  [5],  neglecting  of  the  dissociation  of  nitrogen  should  not  lead  to  a 
significant  error  and  is  reasonable  for  our  simple  model. 

Since  e  (p  is  usually  much  larger  than  the  kinetic  energy  of  an  ambient  gas  molecule  at  the  entrance  plane  of  the 
Hall  thruster  (for  nitrogen  molecules  it  is  about  8.75  eV),  the  npt  can  be  written  as 


nPi  =ne  = 


V0  'Vv 


2  e  <p 

V  M 


(7) 


Substituting  Eq.  (7)  and  V0  as  7.8  km/s  into  an  Eq.  (5),  we  obtain  that  condition  r0  »  Tionlz  is  fulfilled  when 


47.6-  Lnsasaloriiz  (Je ) 
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where  Te  in  Eq.  (8)  is  measured  in  eV. 
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On  the  other  hand,  the  electron  lifetime  in  the  Hall  thruster,  Te_tije ,  must  be  larger  than  the  ionization  time  for  an 
electron,  (VTe  -ngas  <7ipm- )  1 ;  otherw  ise,  the  electron  will  leave  the  thruster  chamber  before  producing  any  ionization 
events;  te  ^  can  be  estimated  as 


L  L 

T"m~V^rip  ~  Ml  E 

It  worth  noting  that  in  Eq.  (9)  we  have  assumed  the  classical  electron  transport  across  the  magnetic  field,  and 
neglected  the  collisions  of  electrons  with  neutral  particles  (since  a  fully  ionized  plasma  is  assumed);  recall  that  in  the 
case  of  a  high  magnetic  field,  electron  transport  in  Hall  thrusters  can  approach  the  classical  limit  [6],  Substituting  in 
Eq.  (9) 
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Me= — - — ,  vej  =2.1- 10-12  -  In  A  •  n  •  T7itZ ,  (0,Ie  =1.76  107  B,  E  =  <p/L,  (10) 

me  ■  vel 


and  taking  into  account  that  (oytelvei  in  Hall  thrusters  is  much  larger  than  one,  we  obtain 


^ e-life 


=  1.5  1015 


_ If  B2 _ 

■J<p- A  •  In  A  •  ngas  -  7^3/ 2 
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where  B  is  measured  in  Gauss,  and  Te  in  eV.  Now,  taking  into  account  that  Te_^e ,  Eq.  (11).  must  be  larger  than 
(VTe  •  nKUS  ■  <7iprU- )”'  we  obtain  that  the  following  condition  on  L  must  also  be  satisfied  to  fulfill  assumption  2; 


41010Ltf-7; 
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Thus,  we  have  shown  that  if  the  conditions  given  by  Eqs.  (8)  and  (12)  are  well  satisfied,  then  the  gas  in  the  air 
breathing  Hall  thruster  is  well  ionized. 

Now  let  us  examine  the  third  assumption  of  the  model-  that  the  plasma  is  quasineutral.  Substituting  np,  from 
Eq.  (7)  into  an  equation  for  the  Debye  radius  [7], 


rD  =  7.43  103 


(  V/2 

T 

\  Pl  7 


(13) 


and  assuming  that  the  gas  is  fully  ionized,  we  obtain  the  following  equation  for  the  Debye  radius: 

V/2 


ro  ~  !04  - 


Tef[<P 

yngas-JA 


(14) 


Thus,  if  the  Debye  radius,  Eq.  (14),  is  much  smaller  than  the  characteristic  dimension  of  the  thruster, 

V/2 
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the  plasma  in  the  Hall  thruster  can  be  considered  quasineutral;  Te  in  Eqs.  (13)  -  (15)  is  measured  in  electron  volts. 

Assumption  4  of  the  model  requires  that  the  characteristic  electron  cyclotron  radius  be  much  smaller  than  the 
thruster  gap.  Substituting  in  an  expression  for  electron  cyclotron  radius  Vre  for  electron  transfer  velocity,  we  obtain 


rb 


-ru  _v  ( rb-ra)-B 


(16) 


where  the  electron  temperature  is  in  electron  volts  and  the  magnetic  field  is  in  Gauss. 

Assuming  that  all  assumptions  of  the  model  are  fulfilled,  the  thrust  of  an  air  breathing  Hall  thruster  can  be 
estimated  as 


Thrust  =  7I-(r£  -r~)  M  V0 


and  the  electric  power  supplied  to  the  discharge  as 

p  =  n  ■  ( rb2  -  r2 )  •  {jr  +  j; )  ■  <p  =  K  ■  ( rfr2  -  r2 )  •  (e ■  np,  ■  Vrdrifi  +  e  ngas  V0 )  <p  , 


(17) 


(18) 


where  n  (r^-r2)  is  the  opening  area  of  the  thruster,  and  M  -V0- ngas  and  AV  ~^2-e  cp! M  in  Eq.  (17)  are 

correspondingly  the  mass  flux  entering  the  thruster  and  the  increase  in  the  velocity  of  heavy  particles  due  to  their 
acceleration  in  the  Hall  thruster  chamber.  Since  the  third  assumption  of  the  model  is  that  the  plasma  is  quasineutral, 
in  Eq.  (18)  in  the  expression  for  we  have  used  npt  instead  of  ne .  Substituting  in  Eq.  (17) 

K.drif,  -  Me  ■  E  ■  vei 1  c%'  >  npi  fr°m  Eq.  (7),  and  V0=  7.8  km/s  and  using  Eqs.  (10),  we  obtain  that 


-.1/2 


lU/V  i,  •  n- np  ls 

jf  =  6.01  10-35 - - *—  , 

B2  L 

(19) 

y,=1.25  10-,5-nsaJ  , 

(20) 

where  is  in  eV,  B  in  Gauss. 

Assuming  diffuse  reflection  of  gas  molecules  from  the  satellite  walls  [3],  the  drag  force  can  be  estimated  as, 

F<trag  -M  -Vq  ■  nsas  ■  S  .  (21) 

It  should  be  stressed  that  since  the  gas  entering  the  thruster  does  not  produce  drag,  the  first  assumption  of  the  model, 
the  opening  area  of  the  thruster,  n  (rp  -  rj ) ,  is  not  included  in  S  . 

One  of  the  important  parameters  governing  the  discharge  in  a  Hall  thruster  is  the  electron  temperature.  Now  let 
us  obtain  an  equation  for  electron  temperature  for  an  air  breathing  Hall  thruster.  An  energy  conservation  equation  in 
the  Hall  thruster  chamber  can  be  written  as 

X  (rb  ~ru  )  "  if  ‘  npl  '  V e.dri/t  +  e '  ngas  'Vq)'*P~  ^wati  +  ^  jet  +  &  anode  >  (22) 

where  the  left-hand  side  of  Eq.  (22)  is  the  electric  power  supplied  to  the  discharge,  Eq.  (18),  while  the  right-hand  of 
this  equation  describes  the  heat  losses  at  the  Hall  thruster  walls  and  at  the  anode,  and  the  energy  flow  leaving  the 
Hall  thruster  with  the  plasma  jet.  Assuming  no  secondary  electron  emission  from  the  wall,  y  =  0  ,  this  corresponds 

to  the  case  of  minimal  heat  loss  to  the  wall;  equations  for  El°2  ,  Eje, ,  and  El°“cIe  can  be  written  as 
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(23) 

(24) 

(25) 


In  Eq.  (23)  we  have  used  the  Bohm  plasma  sheath  theory  [8]  and  in  Eq.  (25)  have  taken  into  account  the  heat  flux 
that  the  electrons  bring  to  the  anode,  and  neglected  the  (typically  small)  ion  flux  to  the  anode.  The  first  term  in  the 
square  brackets  on  the  right-hand  side  of  Eq.  (23)  describes  the  heat  flux  to  the  wall  due  to  the  recombination 
process  plus  the  kinetic  energy  flux  that  ions  bring  to  the  wall,  and  the  second  term  describes  the  heat  flux  that 
electrons  bring  to  the  wall.  The  plasma  sheath  potential  drop,  ,  and  the  Bohm  velocity,  Vg  8,  are 


2  •  n  ■  me  ^ 

M 


(26) 


(27) 


It  should  be  stressed  that  in  our  model  we  have  actually  divided  ions  into  two  parts:  the  fast  ions  in  the  core  of 
the  thruster  chamber  (these  ions  are  accelerated  by  applied  electric  field  and  do  not  reach  the  wall)  and  the  slow  ions 
at  the  vicinity  of  the  wall  which  reach  the  wall.  The  ions  created  at  the  wall  are  accelerated  by  the  sheath  potential 
and  reach  the  wall,  where  they  recombine  with  electrons  and  then  move  back  into  the  chamber  as  neutrals,  and  then, 
they  are  again  ionized  in  the  vicinity  of  the  wall  and  move  back  to  the  wall.  Since  in  our  model  the  neutral  particle 
ionization  length  is  assumed  to  be  much  smaller  than  L .  these  newborn  ions  do  not  gain  much  energy  from  the 
applied  electrical  field  in  this  cycle.  Therefore,  we  have  neglected  the  kinetic  energy  of  these  ions  at  the  sheath  edge 
as  it  has  been  assumed  in  [8], 

Since  the  ion  and  electron  temperatures  in  the  plasma  jet  (in  the  plume.  Figure  1)  are  typically  a  few  electron 
volts  and  much  smaller  than  the  ionization  cost,  which  is  20  -  40  eV,  we  will  further  neglect  Tijet  and  Te  jel  in  Eq. 

(24).  Substituting  the  reduced  equations  for  E  °Zfr ,  Ejel ,  E^de  into  Eq.  (22),  after  some  algebra  we  obtain  a  final 
equation  for  the  electron  temperature: 
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(28) 


Thus,  knowing  the  parameters  of  the  air  breathing  Hall  thruster  ( L ,  ra,  rb ,  B  ,  <p)  and  parameters  of  the 

ambient  air  ( A  ,  criom.  (Te ) ,  agus ,  ngas ,  Tglis ,  E^t )  using  Eq.  (28)  we  can  calculate  Te ,  then  check  the 

assumptions  made  in  the  model.  Eqs,  (1),  (8),  (12).  (15)  and  (16),  and  finally  calculate  the  thrust,  Eq.  (17),  the 
power,  Eq.  (18),  the  ion  and  electron  currents,  Eqs.  (19)  and  (20),  and  other  parameters  of  the  air  breathing  Hall 
Effect  thruster  at  a  given  satellite  orbit 


III.  Numerical  Results 

In  the  numerical  results  presented  in  this  section  the  collision  cross  section  for  molecular  nitrogen  gas.  agas ,  was 
taken  as  4.410  19  m2  [9];  this  corresponds  to 
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The  ionization  cross  section  of  molecular  nitrogen  as  a  function  of  electron  temperature,  Figure  2,  was  calculated 
using  data  [5],  and  In  A  was  taken  as  13.5;  the  length  and  the  inner  and  outer  radii  of  the  Hall  thruster  chamber  were 
chosen  as  0.5,  0.03,  and  0.05  m  respectively.  The  calculations  have  been  performed  for  two  applied  voltages,  3  and 
30  kV,  and  for  satellite  orbits  of  80  and  90  km. 


Te[eV] 

Figure  2.  Ionization  cross  section  of  molecular  nitrogen 

Figure  3  shows  the  drag  factor  fdrag  =  TvaU  V0/  L  calculated  for  80,  90,  100,  and  1 10  km  satellite  orbits.  As  one 
can  see,  the  first  assumption  of  the  model  is  well  fulfilled  only  for  the  80  km  satellite  orbit  However,  since  the 
thrust  of  air  breathing  Hall  thrusters,  Eq.  (17),  is  generally  larger  than  the  maximal  drag  n  ■  (r62  -r*)M  ngas  V02 

that  might  be  produced  by  the  air  flowing  through  the  Hall  chamber,  this  assumption  is  not  critical  for  the  air 
breathing  Hall  thruster  concept.  It  is  worth  noting  that  the  gas  flow  through  the  thruster  is  free-molecular,  Eq.  (2), 
for  H  =  90  km  and  larger,  and  is  collisional.  Eq.  (3),  for  H  =  80  km  and  smaller.  Obviously,  for  different  thruster 
chamber  length  and  width  the  transition  between  the  free  molecular  and  the  collisional  flows  can  take  place  at 
different  orbits. 


80  90  100  110 

H  [kml 

Figure  3.  Drag  factor  for  the  air  breathing  Hall  thruster: 
L  -  0.5  m,  rb  =  0.05  m,  and  ra  =  0.03  m 
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Figures  4  and  5  show  the  ratios  of  r0  / rioniz ,  re_(,yf  VTe  ■  ngas  -  ajuniz ,  and  (rh  -  ra)/rL ,  Eqs.  (8),  (12),  and  (16), 

calculated  for  the  satellite  orbits  of  90  and  80  km  respectively,  the  ionization  cost  of  40  eV.  and  the  applied  voltages 
of  3  kV  and  30  kV;  they  correspond  to  the  second  and  forth  assumptions  of  the  model.  As  one  can  see  for  the 
satellite  orbit  of  90  km  these  assumptions  are  satisfied  for  magnetic  fields  in  the  ranges  of  10  -  30  and  20  -  70  Gauss 
correspondingly  for  <p  =  3  kV  and  30  kV;  and  for  the  orbit  of  80  km  in  the  ranges  of  30  -  120  G  for  <p  =  3  kV  and 
for  (p=  30  kV  in  all  region  of  magnetic  fields,  50  -  300  G.  The  third  assumption  of  the  model,  Eq.  (15),  and 
0)He  /  vei  »  1  condition  are  very  well  satisfied  in  all  ranges  of  magnetic  field,  potential,  and  orbit  considered  in  the 
paper. 


5  10  15  20  25  30  35  40  10  30  50  70  90  110 

B  [G]  B  [G] 

Figure  4.  The  r0  I  Tioni.  -  blue  line,  Te_llfe  ■  VTe  ngM.  ahniz  -  green  line,  {r„-ra)irL-  red 
line  calculated  for  H  =90  km  and  <p  =  3kV  -  (a)  and  p  =  30kV  -  (b);  the  broken  line 
correspond  to  unity. 


B  [G]  B  [G] 

Figure  5.  The  r0  lrioniz  -  blue  line,  re.life-Vre  ngas  c Tioni.  -  green  line,  (rb  -  rj/r,  -  red 
line  calculated  for  H  -  80  km  and  q>  =  3  kV  -  (a)  and  <p  =  30  kV  -  (b). 

Since  the  air  breathing  Hall  thruster  plasma  in  our  model  is  assumed  to  be  fully  ionized,  the  thrust,  Eq.  (20),  is 
determined  only  by  the  ambient  gas  number  density  and  the  applied  voltage,  and  the  ion  current,  Eq.  (17),  is 
determined  by  the  ambient  gas  number  density  only,  The  model  gives  ji  =  8.9- 104  and  4.8- 105  AAn‘  for  H  =  90 
and  80  km  respectively;  and  Thrust  =  18.63,  58.94,  100.5,  and  317.9  N  correspondingly  for  //  =  90  km,  (p=  3  and 
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30  kV  and  H  =  80  km,  cp-  3  and  30  kV.  The  obtained  results  are  very  understandable:  with  a  decrease  in  the 
satellite  orbit,  the  ambient  gas  density  increases,  leading  to  an  increase  in  the  ion  current  and  the  thrust;  and  an 
increase  in  the  applied  voltage  leads  to  larger  AV  and,  therefore,  larger  thrust.  Of  course,  a  lower  orbit  will  lead  to  a 
larger  drag  force  that  the  thrust  will  need  to  compensate  for.  The  drag  will  depend  on  the  satellite  geometry;  this 
paper  considers  the  generic  satellite,  thus,  specific  drag  estimations  go  beyond  the  scope  of  this  paper. 
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Figure  6.  Parameters  of  the  air  breathing  thruster  for  H  -  90  km:  plots  (a),  (e)  and 
(e)  correspond  to  <p  =  3kV  and  plots  (b),  (d)  and  (f)  to  <p  =  30  kV  . 
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Figures  6  and  7  show  the  electron  temperature,  power,  and  _/,•  /  je  vs.  magnetic  field  for  H  =  90  and  80  km 
respectively.  In  these  figures  we  have  selected  magnetic  fields  for  which  the  model  assumptions  are  satisfied.  As 
was  expected  with  an  increase  in  the  applied  voltage,  the  electron  temperature  increases,  and  with  an  increase  in  the 
magnetic  field  it  decreases,  as  shown  in  Figures  6a,  6b  and  7a  and  7b.  Since  the  ion  current  in  the  model  is 
independent  of  the  magnetic  field  (ions  are  not  magnetized),  while  the  electron  current  decreases  with  an  increase  in 
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Figure  7.  Parameters  of  the  air  breathing  thruster  for  H  —  80  km:  plots  (a),  (c)  and 
(e)  correspond  to  <p  =  3kV  and  plots  (b),  (d)  and  (f)  to  <p  =  30  kV  . 
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the  magnetic  field  (they  are  magnetized),  the  electrical  power  decreases  and  the  ratio  of  jL  to  je  increases  with  an 
increase  in  B  ,  as  illustrated  in  Figures  6c  -  6f  and  7c  -  7f.  Since  the  electron  component  of  the  total  current  sharply 
decreases  with  magnetic  field,  as  shown  in  Figures  6e,  6f,  7c,  7f,  P  flattens  for  large  magnetic  fields,  Figures  6c, 
6d,  7c,  7d. 

We  have  also  examined  the  air  breathing  Hall  thruster  with  the  selected  L,  ru ,  and  rh  for  the  satellite  orbits  of 
100  and  1 10  km.  Although  we  could  not  select  B  and  <p  to  satisfy  the  model  assumptions  at  such  high  altitude,  we 

believe  that  changing  the  geometry  of  the  thruster,  for  example  increasing  the  length  of  the  thruster,  see  Eqs.  (8)  and 
(12),  may  allow  this. 


IV.  Concluding  Remarks 

The  idea  of  using  air  breathing  electrical  propulsion  for  a  vehicle  flying  at  orbital  speed  on  the  edge  of  Earth's 
atmosphere  has  been  examined  for  a  thruster  based  on  the  Hall  Effect.  We  have  shown  that,  conceptually,  such  a 
thruster  can  indeed  work  effectively  at  the  orbits  of  about  80  -  90  km,  producing  significant  thrust  in  the  range  of  20 
to  320  N  for  considered  conditions. 

It  should  be  pointed  out  that  the  power  required  maintaining  such  high  level  of  thrust  is  in  the  order  of  MW  with 
thrust-to-power  ratio  of  4.37  -  13  |iN/kW  which  is  comparable  with  a  typical  high-efficient  electric  propulsion 
device.  Recall  that  while  air-breathing  concept  seems  to  be  reasonably  achieved  in  terms  of  thrust  production  it  can 
be  realized  with  high  power  available  only.  Thus,  analysis  of  available  power  generation,  power  transmission  is 
waiTanted  as  a  part  of  the  overall  air-breathing  thruster  and  satellite  analysis. 

Note,  that  this  paper  has  considered  relatively  simplified  model  of  the  plasma-wall  interaction  inside  the  thruster 
chamber.  In  particular,  zero  secondary  electron  emission  (SEE)  scenario  was  invoked  to  describe  the  sheath 
formation.  Such  treatment  is  limited  to  the  wall  materials  with  a  low  SEE  [10].  Further  analysis  of  the  air-breathing 
concept  should  involve  broader  conditions  for  the  SEE  thus  extending  operational  parameters  of  this  device. 

We  believe  that  by  changing  the  thruster  geometry,  we  could  find  reasonable  conditions  allowing  increasing  the 
flying  orbit  of  an  air  breathing  satellite  up  to  100  -  1 10  km.  It  should  be  stressed  that  in  the  model  we  have  assumed 
that  the  plume  is  fully  ionized  corresponding  to  the  maximum  possible  thrust  and  power.  However,  if  we  soften  this 
condition,  assuming  that  the  plume  can  be  only  partially  ionized,  the  parameters  of  the  thruster  become  less 
“extreme"  and  may  lead  to  use  of  the  thruster  at  higher  orbital  altitude. 
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